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ABSTRACT 



Context. The next generation near- and mid-infrared Galactic surveys will yield a large number of new highly obscured star clusters. 
Detailed characterization of these new objects with spectroscopy is time-consuming. 

Aims. Diagnostic tools that will be able to characterize clusters based only on the available photometry will be needed to study large 
samples of the newly found objects. 

Methods. The brightness difference between the red clump and the main-sequence turn-off point have been used as a model- 
independent age calibrator for clusters with ages from a few 10** to 10'" yr in the optical. Here we apply for the first time the method 
in the near-infrared. 

Results. We calibrated this difference in A'-band, which is likely to be available for obscured clusters, and we apply it to a number of 
test clusters with photometry comparable to the one that will be yielded by the current or near-future surveys. 

Conclusions. The new calibration yields reliable ages over the range of ages for which the red clump is present in clusters. The slope 
of the relation is smoother than that of the corresponding V-band relation, reducing the uncertainty in the age determinations with 
respect to the optical ones. 

Key words. Open clusters and associations : general - Milky Way : structure - Infra Red surveys. 



1. Introduction 

The recent advance in the infrared (IR) instrumentation made it 
possible to carry out large scale sky surveys. The Two Micron 
All Sky Survey (2MASS; Skrutskie et al. 2006) was the first one 
to provide a deep look into the most obscured regions of the 
Milky Way. A number of new clusters were found in 2MASS 
(i.e. Dutra & B ica|2000 Ivanov et al. |2002bl Bica et al. |2003l 
Borissova et al. 12003] ). Currently, th e UKI RT Infrared Deep Sky 
Survey (UKIDSS; Lawrence et al. I2007I I is creating a map of 
selected Northern sky regions, typically ~5 mag deeper than the 
2MASS, and in the near future the Visible and Infrared Survey 
Telescope for Astronomy (VISTA; Emerson et al. 120061 Dalton 
et al. 2006) will provide similar coverage in the South. These 
projects include specialized multicolor Galactic plane surveys: 
UKIDSS Galactic plane surveys (GPS; Lucas et al. |2008l) and 
VISTA variables in the Via Lactea (VVV; Minniti et al. |2006] ) 
which are expected to reveal a large number of new star clusters. 

Their characterization is a time-consuming process, requir- 
ing spectroscopy in the near-IR, because of the strong inter- 
vening dust absorption (Aa:~0.1xAv; Rieke & Lebofskv ll985l ). 
This prompted us to attempt to develop various methods to deter- 
mine the cluster parameters from the photometry alone. For ex- 
ample, the metal abundances of globular clusters is easily mea- 
sured with an error of cr([Fe/H])~0.2 dex (i.e. Ferraro et al. |19991 
Ivanov & Borissova '2002a') from the slope of the RGB on the 
color-magnitude diagrams (CMD). 

In this Note we focus on a new age calibrator based on the K- 
band. difference between the clump and the turn off point (TO) 



in the CMD of star clusters older than about 300 Myr, which 
will be widely used to get a quick and approximate age values 
for newly discovered IR star clusters. 

The layout of the paper is as follows. In Section 2 we in- 
troduce the concept of age calibrator as derived and used in the 
optical regime. Section 3 describes the sample of very well know 
star clusters on which we base our new age calibrator. A few an- 
alytical relations are given in Section 4. These relations and then 
applied to a few test-cases in Section 5. Section. 6, finally, sum- 
marizes out findings. 



2. Age calibrations for open star clusters 

The magnitude difference AV between the main sequence (MS) 
TO and the clump of He burning stars is a well known age indi- 
cator for Galactic open clusters. The very requirement to have a 
well-defined red clump limits it to clusters older than ~300Myr 
when the He burning stars appear. However, the method is dis- 
tance and reddening independent, as demonstrated early on by 
Cannon ( 19701 ), and it is not surprising that it has been widely 
used ever since (see, for instance, Carraro & Chiosi i 19941 Salaris 
et al. ,2004j . 

Different, more sophisticated versions of this caUbrator have 
been devised in the past, in the optical regime, by Anthony- 
Twarog & Twarog ( 119851 ) and Phelps et al. (1994). They also 
include the color difference between the TO and the red giant 
branch (RGB). However, their application is more limited be- 
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Table 1. Parameters of the clusters' sample. The last column 
gives a reference to the literature source of the adopted cluster 
age and metalhcity. 



M 67 



NGC 2506 



Name 


RA 


DEC 


Age 


\Feim 


aa: 


Ref. 




hh:mm:ss 




[Gyr] 




[mag] 




NGC 188 


00:47:28 


+22.384 


6.0 


+0.01+0.08 


4.0+0.15 


1 


NGC 752 


01:57:41 


-23.254 


1.6 


-0.16±0.10 


2.7+0.20 


2 


NGC 2420 


07:38:23 


+ 19.634 


2.1 


-0.26±0.07 


3.3+0.20 


3 


NGC 2477 


07:52:10 


-05.850 


0.8 


+0.07±0.03 


1.9+0.25 


4 


NGC 2506 


08:00:01 


+09.935 


1.9 


-0.37+0.02 


3.1+0.25 


5 


NGC 2509 


08:00:48 


+05.820 


1.2 


+0.00+0.10 


2.5+0.20 


6 


M67 


08:51:18 


+31.896 


4.0 


+0.00+0.01 


3.8+0.10 


7 


NGC 3680 


11:25:38 


+ 16.918 


1.8 


-0.16+0.03 


2.9+0.15 


8 


NGC 6 134 


16:27:46 


-00.200 


1.0 


+0.18+0.10 


2.0+0.15 


9 


NGC 6253 


16:59:05 


-06.260 


3.0 


+0.36+0.07 


3.7+0.15 


10 


IC4651 


17:24:49 


-07.907 


1.1 


-0.11+0.01 


2.1+0.20 


11 


NGC 6819 


19:41:18 


+08.481 


2.5 


+0.09+0.03 


3.4+0.25 


12 


NGC 6939 


20:31:30 


+ 12.304 


2.2 


+0.02+0.05 


3.2+0.15 


13 


NGC 7789 


23:57:24 


-05.385 


1.4 


-0.10+0.05 


2.4+0.20 


14 



1) Fomal et al. |2007l 2) Anthony-Twa rog & Twarog l2006bl 3) Anthon y- 
Twarog et al. 2006a: 4) Kassis et al. fTWTl 5) Ma rconi et al. fTWTI 6) 
Carraro & Costa 2007 ; 7) Sarajedini et al.l2009l, 8) Antho ny-Tw arog 
& Twai-og 20 04, 9 ) Bruntt et al. 1999: 10) Twarog et al. l2003l 11) 
Meibon et al. i2002[ 12) Kalirai et al. 120011 13) Andreuzzi et al. l2004l 
14) Vallenarietal. l2000l 



cause usually the RGBs of open clusters are poorly populated, 
and the RGB colors are ill-defined. 

The simple difference AV between the red clump and TO 
- the two major features in the CMD - has the advantage of 
providing a quick age estimate even from un-calibrated CMDs 
(Phelps et al. [19941 Salaris et al. |2004] i. As described by Carraro 
& Chiosi (1994) this indicator shows a slight metallicity depen- 
dence in the optical: AV at fixed age is larger for higher metal- 
licity clusters than for lower metallicity ones. 

The aim of the present work is to extend this calibrator to the 
near-IR, where the metallicity effect should be weaker. We will 
use the 2MASS filter as the least affected by the extinction, 
and throughout the Note will mark it as K, for simplicity. 



3. Sample Selection and Indicator Definition 

We selected 14 open clusters with 2MASS photometry, with reli- 
able age and metallicity estimates, spanning a wide range of the 
parametric space (Table[T]). The data come from the latest edition 
of the WEBDAQ database maintained by E. Paunzen at Vienna 
University. Throughout the paper we adopted a conservative es- 
timates of the uncertainty in the age of 15% and a factor of 2 for 
the abundance [Fe/H]. 

Some of the clusters suffer from significant fore- and back- 
ground contamination since they are located low onto the 
Galactic plane. To minimize it, we built up homogeneous CMDs 
considering only the stars within the accepted cluster radii es- 
timates from Dias et al. (i2002| i. The CMDs for M67 and 
NGC 2506 - two clusters with rather different ages and abun- 
dances - are shown in Fig. 1. 

Next, we identified the clump and TO magnitudes, re- 
spectively at /:^'""''^=8.00±0.05, and 7i:^^=11.80±0.10mag, 
for M67. This implies A/:=3.8+0.1 mag. For NGC 2506 we 
found /:'''™P=10.85+0.15, and /i:™=14.00+0.20mag, yielding 
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Fig. 1. Definition of the parameter A/T. The left panel shows 
the CMD of M 67 (solar metallicity, 4 Gyr old), the right shows 
the CMD of NGC 2506 (significantly subsolar metallicity, twice 
younger than M 67). 
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Fig. 2. Age vs. t^K relation for the clusters listed in Table[T]ln 
the same Table we report the uncertainties in age and A/T. 



A/r=3. 15+0.25 mag. These examples illustrate that for an old 
cluster like M 67 the definition of the clump magnitude is easier 
despite the smaller number of He burning stars, while the clump 
at younger ages is wider both in color and magnitude due to evo- 
lutionary effects (Girardi et al. 2000). 

The measured values of KK are also listed in Table[T] Their 
uncertainties were obtained tentatively, by eye. These take 
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Fig. 3. Left panel: Age vs AK relation for the 14 clusters listed in 
Table[T] The polynomial fit (Eqn.|2]i is showed as a dashed line, 
and the linear one (Eqn.[T]l - as a solid line. Right panel: The 
optical relation from Carraro & Chiosi (1994). The sohd line is 
a linear fit through the data. 

into account the broadening of features in the CMD due to 
photometric errors, patchy extinction and unresolved binaries. 
The experience shows that such estimates can be even more 
rehable than the automated ones (i.e. see the discussion in 
Section. 4 of Phelps & Janes 1994). The age (logarithm) versus 
AK relation is plotted in Fig.|2] 

We emphasize that this IR relation looks smoother than the 
optical counterpart (Carraro & Chiosi 1994) mainly because 
of the significantly weaker metallicity dependence (see below). 
Furthermore, the smooth points distribution in Fig. 2 makes any 
analytical fit through the data easier and more robust . 

4. Calibration 

An analytical representation of the age versus AA" relation is 
needed to facilitate the derivation of cluster ages. The relation 
appears non-linear, prompting us to fit a higher order polyno- 
mial: 

log(Age) = 9.09(0.29) + 0.24(0.20)xA/:-h0.10(0.03)x(AA:)2(1) 

with r.m.s.=0.28. The reduced of the fit is 0.08. 
For comparison, the linear fit to the data is: 

logiAge) = 8.27(0.1) + 0.35(0.03) x AK (2) 

with r.m.s.=0.36. Numbers in brackets are corresponding uncer- 
tainties. The reduced of this fit is 0.04. 

The age in these relations is expressed in Gyrs. The two fits 
are nearly indistinguishable for 2.5<A/r<3.5 mag - the devia- 
tions from the linearity occur in fact at the extremes of the age 
range. The fits are plotted in the left panel of Fig. [3] while in the 
right panel we show the optical relation from Carraro & Chiosi 
(1994), for the sake of comparison. 

As anticipated in the previous Section, the data point distri- 
bution in the IR plot (left panel) is smoother than the optical one 
(right panel), and allow a better fit through the points. 

Next, we include the metallicity in the fit, obtaining a 3- 
dimensional relation: 



I 




Fig. 4. Age vs AK for the 14 clusters listed in TablefT] including 
the metallicity dependence. The projection on the age vs AK is 
also shown. 



As expected, the dependence on metallicity is negligible, 
making AK very useful, since in most cases the newly dis- 
covered clusters will not have metallicity estimates available. 
In comparison with the optical relation (Carraro & Chiosi 
1994 Eqn. 2), the linear relation in the IR is less sensitive to 
metallicity, and also shows a somewhat shallower dependence 
on AK (0.34 versus 0.45, see CaiTaro & Chiosi 1994). 

As for the uncertainty on the age, it is easy to show - by simple 
propagation - that in the linear cases (Eqns. 2 and 3), the uncer- 
tainty in the age does not depend on the value of AA'. This is not 
the case of course for Eqn. 2, which has a quadratic dependence 
on AK. 



5. Validation of the new calibrations 

To demonstrate the power of the new calibrations, we applied 
them to estimate the ages of three clusters with shallow 2MASS 
photometry, but that have IR photometry from other sources: 
NGC6791, NGC2141, and Berkeley 17. Spectroscopic metal- 
licity estimates are also available for them, allowing us to use 
all three relations in a metallicity range which encompasses the 
typical metal content of Galactic open clusters. The results are 
listed in Table|2] The age uncertainties were derived propagating 
the measurement errors through the calibrations. 

5.1. NGC6791 

IR photometry for this cluster is taken from Carney et al. ( I2005I I. 
It was acquired with the IRIM camera at the 4-m Mayall tele- 
scope at Kitt Peak National Observatory. The data were pho- 
tometrically calibrated with the UKIRT standards (Casali & 
Hawarden |1992| l, which were also used to calibrate the 2MASS 
itself. Not surprisingly, Eqn. 37 from Carpenter ( 120011) : 

KiMAss = /:f//f/R7-+0.004(0.006)x(7-7i:)t/^/«7-+0.002(0.004)(4) 



logiAge) = 8.28(0.08)H-0.34(0.03)xA/:+0.00005(0.1)x[Fe/i/](3) indicates that the color difference A(7-/:)~0.3 mag between the 

TO and clump leads to only 0. 0012+0. 0018mag difference be- 
with r.m.s.=0.06 (Fig.lUl. The reduced of the fit is 0.004. tween the K-band magnitudes of the two systems - much smaller 
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than our measurement errors. The constant term is neghgible system. We provide Hnear and second order polynomial fits of 
for calculating the AK parameters, and it is statistically indis- the age versus AK relationship, and a multi-linear fit which also 
tinguishable from zero, anyway. takes into account the metallicity of the clusters, when available. 

The CMD of NGC 6791 is shown in the left panel of Fig.|5] However, it is unlikely that the metallicity will be immediately 
Widely accepted values for age and metallicity are in the range of available for most of the new objects, so Eqns. 1 and 2 will be 
7-9 Gyrs and [Fe/H]=-i-0.35 to +0.45 dex, respectively(Carney et mostly used used for en-masse analysis of the new clusters, 
al. 2005, and references therein). Our equations 1, 2, and 3, yield Beside, we show that the dependence on metallicity is very 
ages of 9.3, 6.5 and 6.0 Gyrs, respectively. shallow, which makes these relations of more robust use than 

their optical counterpart. 

5 2 NGC 2141 ^^^^ illustrate with three examples that the simple linear 

■ ■ relation works better for intermediate AK values, while the 

IR photometry for this cluster is taken from Carraro et al. (2O0T). polynomial relation is preferable at the exti-emes of the validity 
It was acquired with the ARNICA camera at the Tirgo telescope, range. 
The data were photometrically calibrated with the ARNICA 

standards (Hunt et al. il998i) . Eqn. 8 from Carpenter COOlt : The method can potentially be useful for a wider age range, and 

these calibrations should be extended in the future to cover both 
KiMASs = /Taka'/ca -0.024(0.01 l)x(7-/:)AfiA'/cA+0.012(0.006)(5) younger and older clusters. 



indicates that the color difference A( J — K)~0.25 mag between 
the TO and clump leads to only ~ 0.006+0.003 mag difference 
between the K-band magnitudes of the two systems - again, 
much smaller than our measurement errors. 

The CMD of NGC 2141 is shown in the middle panel of 
Fig.|5] For this cluster, metallicity is [Fe/H]=-0.26, and the age 
is ~2.5Gyrs (Carraro et al. 2001). We obtained ages of 2.2, 2.4 
and 2.3 Gyrs, respectively from Eqns. 1, 2, and 3. These re- 
sults are in good agreement with the isochrone fitting estimates 
(Carraro et al. lIOOTT i. 

5.3. Berkeley 1 7 

IR photometry for this cluster is taken from Carraro et al. (|1999l l. 
It was acquired with the ARNICA camera at the Tirgo tele- 
scope, and photometrically calibrated with the ARNICA stan- 
dards (Hunt et al. 1998). As shown in the previous subsection, 
our age measuring parameter AK is not affected significantly by 
using a different photometric system than the 2MASS. 

The CMD of Berkeley 17 is shown in the right panel of 
Fig.m This cluster is known to have marginally subsolar abun- 
dance [Fe/H]=-0.10 dex and is it 9-lOGyrs old (Friel et al. 
2005). Our calibrations yielded ages of 10.0, 7.6 and 7.0 Gyrs, 
respectively for Eqns. 1, 2, and 3. 

This exercise shows that the age of a cluster within the lin- 
ear regime of our calibrator (2.5<A^r<3.5 mag) is closer to the 
isochrone-derived one when simple linear (Eqn. 1), or the metal 
dependent (Eqn. 3) fitting are used. The most deviating relation 
is the polynomial one. In the other two cases of very old clusters, 
the polynomial relation (Eqn. 2) is the most accurate, and pro- 
vides age estimates very close to the widely accepted in the liter- 
ature. Finally, the new calibrations provide, in the linear regime, 
age estimates that agree quite well with isochrone-based ages, 
while outside this range the difference might be of the order of 
30%. 

6. Summary and conclusions 

In this paper we described relations of the magnitude difference 
between the TO and the clump AK in the IR CMD versus age for 
800-1000 Myr old open star clusters. The age limits are set by the 
life time of the red clump stars and available photometric data. 
Our goal is to develop a simple tool to measure the age of newly 
detected star clusters from the next generation IR sky surveys. 
The relations are derived in the widely used 2MASS photometric 
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Fig. 5. Application of the new age indicator to NGC 6791, NGC 2141 and Berkeley 17 (from left to right). AK is marked on each 
CMD. 

Table 2. New age estimates for the three test-cases discussed in this paper. The last three columns list the ages determined from 
Eqns. 1, 2 and 3, respectively. References for the literature ages and abundances are given in the text. 



Name RA 


DEC 


Age 


[Fe/H] 


AK Age{l) 


Age(2) 


Age(3) 


hh:mm:ss 




[Gyr] 




[mag] [Gyr] 


[Gyr] 


[Gyr] 


NGC 6791 19:20:53 


+37:46:18 


8.0 


+0.40 


4.40±0.20 9.3±1.3 


6.5+0.5 


6.0+0.5 


NGC 2141 06:02:55 


+10:26:48 


2.5 


-0.26 


3.20+0.20 2.2+0.2 


2.4+0.2 


2.3+0.2 


Berkeley 17 05:20:36 


+30:36:00 


9.0 


-0.10 


4.60+0.20 10.0+1.4 


7.6+0.6 


7.0+0.7 



